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MINIMUM-MASS DESIGN OF FILAMENTARY COMPOSITE PANELS 
UNDER COMBINED LOADS: DESIGN PROCEDURE BASED 

ON A RIGOROUS BUCKLING ANALYSIS 

W. Jefferson Stroud, Nancy Agranoff, and Melvin S. Anderson 
Langley Research Center 

SUMMARY 

A procedure is presented for designing uniaxially stiffened panels made of composite 
material and subjected to combined inplane loads. The procedure uses a rigorous buckling 
analysis and nonlinear mathematical programing techniques. 

Design studies carried out with the procedure consider hat-stiffened and corrugated 
panels made of graphite -epoxy material. Combined longitudinal compression and shear 
and combined longitudinal and transverse compression are the loadings used in the studies. 
The capability to tailor the buckling response of a panel (i.e., design a panel so that it will 
have specified buckling loads at various buckling wavelengths) is also explored. Finally, 
the adequacy of another, simpler, analysis-design procedure is examined. 

The report demonstrates that a panel design procedure with a high-quality buckling 
analysis and with complete generality of constraints is practical. Such a procedure can 
be used to avoid failure from complex buckling modes and to determine mass and propor- 
tions of panels for multiple design load conditions and constraints. 

INTRODUCTION 

Stiffened panels made of metal and/or composite materials have wide application in 
aerospace structures. In an effort to increase the structural efficiency of these panels 
and, at the same time, to account for numerous manufacturing and other design require- 
ments, panels incorporating advanced design concepts are being explored. (See, for 
example, refs. 1 and 2.) Because these panels exhibit complex buckling modes, it is often 
necessary to use relatively sophisticated analyses to examine their stability. In recent 
years, concurrent advances in analysis techniques (see, for example, refs. 3 to 6) and 
synthesis concepts (see, for example, refs. 7 to 10) have made it possible to develop 
analysis-design procedures that can cope with the complex buckling modes and, thereby, 
exploit the potential of the advanced structural design concepts. One such analysis-design 
procedure is described in this report. 



This report presents a procedure for designing uniaxially stiffened panels having an 
arbitrary cross section, made of composite material, and subjected to inplane longitud- 
inal N x , transverse Ny, and shear N X y loadings. The procedure is based on a rigor- 
ous buckling analysis and nonlinear mathematical programing techniques. Design studies 
which consider graphite -epoxy hat- stiffened panels and graphite -epoxy corrugated panels 
are also presented. In these studies, mass-strength charts are developed for panels sub- 
jected to combined longitudinal compression and shear and combined longitudinal and 
transverse compression. The capability to "tailor" the buckling response of a panel (i.e., 
to design a panel so that it will have specified buckling loads at various buckling wave- 
lengths) is also explored. Finally, the adequacy of a simpler analysis -design procedure 
(ref. 11) is examined. 

SYMBOLS 

Values are given in both SI and U.S. Customary Units. The calculations were made 
in U.S. Customary Units. 

A surface area of panel, pL 

b^ element lengths defined in figure 8 

C constant defined in equation (Bl) 

Di6 ? D 2 6 laminate stiffnesses relating bending behavior with torsional behavior 
laminate bending and twisting stiffness matrix 
d distance between stiffeners, 2b^_ 

E 1; E 2 Young's modulus of composite material in fiber direction and transverse to 
fiber direction, respectively 

EZ1,EZ2,EZ3 eccentricities or offsets in the z-direction shown in figure 19 

F^y^ buckling mode shape of panel element in y^-direction 

Gj 2 shear stiffness of composite material in coordinate system defined by 

fiber direction 
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integers 


hi 

K constant defined in equation (B2) 

k element number 

L panel length in x-direction (see fig. 1) 

required local buckling load for longitudinal compression loading 

Nq required overall buckling load for longitudinal compression loading 

N x applied longitudinal compressive loading per unit width of panel (see fig. 1) 

N value of N that causes buckling; may be for specified X 

x cr x 

^N x ^ value of (N x ^ that causes local buckling of element k 

(^x)design va -l u e of N x for which panel is designed 

(N x )k longitudinal compressive loading in element k 

N X y applied shear loading per unit width of panel (see fig. 1) 

N value of shear load that causes buckling 

xy cr 

/N„. \ value of ( N \ . that causes local buckling of element k 
\ ^cr/k \ xy / K 

k shear loading in element k 

Ny applied transverse (y-direction) loading per unit width of panel 

n integer 

p period of stiffened panel, b^ + 2b^ 

Re( ) real part of complex quantity within parentheses 
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t 2 or t 2 (±45) thickness of ±45° sheet 

t^(0) thickness of lamina in element k at fiber orientation angle 0 (see fig. 8) 

u,v,w buckling displacements (see fig. 17) 

W mass of one period of stiffened panel 


mass index 


X,Y,Z axes (defined in figs. 1, 17, and 18) 


x,y,z 


coordinates in longitudinal, transverse, and lateral directions, respectively 


local coordinates for each element making up a panel 


design variables; values of the design variables for an arbitrary design 
values of the design variables at the initial point of a Taylor series expansion 

a a a 

e V e 2’ y l2 a ll° wab le lamina strains in material coordinate system defined by fiber 
direction 

6 fiber orientation angle shown in figure 1 

X buckling half -wavelength, L/n 

M 12^21 Poissons ratios of composite material in coordinate system defined by fiber 
direction 


density 


a a a 

a l ,cr 2’ T 12 a U° wa ble lamina stresses in material coordinate system defined by fiber 
direction 
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Subscripts: 


cr value that causes buckling 

k element number 

i local 

A subscript preceded by a comma denotes partial differentiation with respect to 
the subscript. 


ANALYSIS -DESIGN PROCEDURE 

The procedure described herein allows a high-quality, eigenvalue buckling analysis 
to be used to design minimum-mass, uniaxially stiffened panels having an arbitrary cross 
section and subjected to combined inplane loads (fig. 1). The procedure consists of three 
basic components: 

(1) A linked-plate, eigenvalue buckling analysis for stiffened composite panels. The 
term "linked-plate" is used in this report to denote the type of buckling analysis used in 
references 3 to 6. 

(2) A nonlinear mathematical programing optimizer for function minimization with 
inequality constraints. 

(3) Approximation concepts such as Taylor series expansions of constraints and 
selective constraint retention. 

These three components and a material strength analysis were combined in a 
research-oriented computer program to study hat-stiffened and corrugated panels. In 
this section of the report, the procedure and, in some instances, the computer program 
that implements the procedure are described. Emphasis is placed on the buckling analy- 
sis and on the approximation concepts associated with the buckling constraints. Initial 
buckling is considered to be a failure mode - that is, it is assumed that the panel has no 
postbuckling strength. The stress analysis and strength criteria used in the computer 
program are similar to those used in a panel design procedure described in reference 11. 

Loadings Considered 

The loadings considered are shown in figure 1. They are: 

(1) A uniform, longitudinal compressive loading N x in the direction of the 
stiffeners. 
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(2) A uniform transverse compressive loading Ny. 

(3) A uniform shear loading N X y. 



Buckling Analysis 

The buckling analysis used in the computer program is an efficient, stiffened-panel 
buckling analysis code denoted VIPASA (Vibration and Instability of Plate Assemblies 
including Shear and Anisotrophy) which is described in references 3 and 4. The VIPASA 
analysis treats an arbitrary assemblage of plates with each plate loaded by N x , Ny, 
and N X y. The response of each plate element making up the stiffened panel is obtained 
using an exact solution of the thin-plate equations. The analysis connects these individual 
plate elements and maintains continuity of the buckle pattern across the intersection of 
neighboring plate elements. An illustration of a continuous buckle pattern is shown in 
figure 2. The stiffened panel is assumed to be uniform in the x-direction (fig. 1) and 
simply supported along the edges x = 0 and x = L. For these reasons the buckle pat- 
terns in the x-direction are taken to be sine waves whose half -wavelengths are fractions 



BUCKLING-MODE SHAPE 

UNDEFORMED SHAPE 

Figure 2.- Continuous buckle pattern. 
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(1/n) of the panel length. In this report, the type of analysis described above is denoted 
a "linked-plate" analysis. A more detailed discussion of VIPASA, including comments 
on the panel modeling used for the hat-stiffened and corrugated panels considered in this 
report, is presented in appendix A. 


Design Procedure 

In the panel design procedure described in reference 11, the buckling loads are cal- 
culated with relatively simple, explicit equations that require only a small amount of 
computer time for evaluation. The analysis is, therefore, coupled directly to the opti- 
mizer. In contrast, in an analysis such as VIPASA the buckling loads are calculated with 
an eigenvalue buckling analysis which is much more time-consuming than the explicit 
equations of reference 11. For that reason, a more sophisticated and efficient sizing 
technique is employed in the present procedure. The most important feature of the sizing 
technique is that it is based on a Taylor series approximate analysis approach (refs. 7 
and 8). The optimizer, denoted CONMIN (ref. 10), used in the computer program is 
based on a feasible direction algorithm. In this approximate analysis approach the opti- 
mizer uses only approximate values of the constraints based on a sequence of Taylor 
series expansions of the constraints. The present strategy greatly reduces the number 
of complete VIPASA analyses required during the synthesis. The general program 
organization is shown in figure 3. 


ANALYSIS 

MODULE 


TAYLOR SERIES 
MODULE 


RESIZING 

MODULE 



Figure 3.- General approach used in present procedure. 
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The Taylor series expansions of the constraints are of the form 


N, 


x 


cr 





(i) 


in which N x is the buckling load, a- are the design variables, and a. are the 

values of the design variables at the initial point of the expansion. The buckling load 
at a j and the derivatives of the buckling load at 5- are evaluated with VIPASA. (For 
the studies presented herein, finite-difference approximations were used to calculate the 
derivatives.) Taylor series expansions are used for the strength constraints as well as 
for the buckling constraints. 

The overall synthesis strategy, which is shown schematically in figure 4, consists 
of a series of subsyntheses in which the optimizer adjusts the values of the design vari- 
ables based on approximate values of the constraints (eq. (1)). An upper limit is imposed 
on the change of each design variable during each subsynthesis to insure the adequacy of 
both the list of constraints that are considered to be active and the Taylor series expan- 
sions of those constraints. These limits to the changes in the design variables are 
referred to as move limits. The move limits associated with each subsynthesis are indi- 
cated by the dashed rectangles in figure 4. (The move limits used in the studies pre- 
sented herein were usually 10 percent of the values of the design variables at the 
beginning of a subsynthesis.) The solid circular symbol at the center of each rectangle 



Figure 4.- Synthesis strategy including move limits during subsynthesis 

in two-design-variable space. 
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represents the point at which the Taylor series expansions are carried out for each sub- 
synthesis. The end point of one subsynthesis becomes the initial point of the next sub- 
synthesis. Accurate values of the constraints and derivatives of the constraints are then 
recalculated, and new Taylor series expansions are generated. Ten subsyntheses are 
usually adequate to obtain convergence if the initial design is reasonably well chosen. 

(For the studies presented herein, typical run times on Langley's CDC CYBER 175 com- 
puter were 100 to 300 seconds for orthotropic panels without shear and 200 to 500 seconds 
for panels that are anisotropic and/or have a shear loading.) 

An example which demonstrates the buckling-constraint aspects of the design prob- 
lems that must be treated with the procedure is presented in figure 5. In this example a 


o CALCULATED POINTS 



0 " .2 ' .4 .6 ‘ .8 1.0 

1 

L 

Figure 5.- Ratio of buckling load to design load as a function of buckling half -wavelength 
for a hat-stiffened panel designed for a loading of N x /l = 689 kPa (100 lbf/in^); 
tg = 0.56 mm (0.022 in.). 

hat- stiffened panel is designed to support a longitudinal compressive loading of 

N x /l = 689 kPa (100 lbf/in2). The transverse and shear loads are taken to be zero. 

The buckling response diagram for the final design is shown in figure 5. In this diagram 

the buckling load is given as a function of the buckling half -wavelength. Both scales are 

nondimensional for convenience. The object of the design procedure is to insure that the 

buckling load N is greater than or equal to the design load ( N 
x cr v 

buckling half -wavelengths X given by A = L, L/2, L/3, L/4, .... In the case shown 

in figure 5, the critical half -wavelengths are A = L, which corresponds to an overall 

buckling mode, and A = L/28, which corresponds to a local buckling mode. The buckling- 

mode shapes for these two buckling modes are given in figures 6 and 7, respectively. In 

this example, buckling modes corresponding to A = L/2 and L/3 are near critical; for 

some cases, these wavelengths became critical at the final design. 
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Figure 6.- Overall-buckling-mode shape for panel whose buckling response diagram 

is shown in figure 5. A = L. 



Figure 7.- Local buckling mode shape for panel whose buckling response diagram 

is shown in figure 5. A = L/28. 

Buckling at each half -wavelength was considered to be a separate constraint. Only 
a few of these buckling constraints need to be retained at any stage of the synthesis. The 
others can be eliminated using various constraint deletion techniques (ref. 8). In the com- 
puter program used to obtain the results presented herein, a buckling half -wavelength A 
was considered to be active if it was prescribed or if it provided the lowest buckling load 
at some stage of the synthesis. 


Cycling or Mode Switching 

There are many buckling eigenvalues at each value of the buckling half -wavelength 
A. All of these buckling eigenvalues can be obtained from VIPASA. Although the 
computer program used to carry out the studies in this report contains sufficient logic 
to account for several simultaneous buckling eigenvalues at their respective buckling wave- 
lengths, the program examines only one - the lowest - of the buckling eigenvalues at a 
given wavelength. A second, higher buckling eigenvalue at that same wavelength may have 
different buckling characteristics and, after a small change in the design, may become the 
lowest buckling eigenvalue at that same wavelength. Unless buckling eigenvalues with both 
characteristics are taken into account simultaneously, the procedure may cycle between 
them and make no progress. Cycling occurred during several of the studies reported 
herein. (Cycling was controlled by treating one or two of the panel dimensions as fixed 
parameters rather than as free variables.) An improvement is needed that would identify 
and provide Taylor series expansions of the critical eigenvalues for more than one buckling 
mode at the same wavelength. 
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DESIGN STUDIES 


Design studies are presented for panels subjected to two types of loadings: 

(1) combined longitudinal compression and shear, and (2) combined longitudinal and trans- 
verse compression. The capability to tailor the buckling response of a panel and the 
implications of this tailoring are also explored. In some cases the results of a simpler 
analysis -design procedure (ref. 11) are given for comparison. Studies which examine the 
adequacy of the simpler procedure are presented in appendix B. The buckling boundary 
conditions along the edges x = 0 and x = L are simple support, and the boundary con- 
ditions along the longitudinal edges (y is constant) are either simple support or sym- 
metric. (The boundary conditions are discussed in more detail in appendix A.) 

Configurations 

The panel configurations considered are the hat-stiffened and corrugated configura- 
tions studied in reference 11 and shown in figures 8(a) and 8(b), respectively. The cross- 
section geometry is defined in terms of the four dimensions bj, b2, bg, and b^. The 
quantities b^, bg, and bg are the widths of elements 1, 2, and 3, respectively. The 
quantity b^ is the half-width of element 4. Because of symmetry, b^ is replaced 
by b 3/ /2 in the corrugated panel. 

The panels are assumed to be made of a filamentary composite material with ply 
orientation angle 6 measured from a line drawn in the direction of the stiffeners (the 
x-direction) as shown in figure 1. In the examples presented in this report, only two 
fiber orientations are considered: 0 = 0° and 0 = ±45°. The 45° plys are assumed to 

be balanced and symmetric and are, therefore, denoted ±45°. For some loadings, in par- 
ticular those involving a transverse load, the addition of 90° plys may provide more effi- 
cient designs; however, 90° plys are not used in the studies presented herein. 

The configurations are defined so that the 0° plys and the ±45° plys have specific 
locations (fig. 8(c)). All elements contain ±45° plys. Elements 3 and 4 can also contain 
0° plys. In both configurations tg(0°) refers to the thickness of one of the two equal 
0° layers in element 3. In the hat-stiffened configuration the thickness of a 0° layer sand- 
wiched between ±45° skins in element 4 is denoted t^(0°). 

The hat-stiffened configuration is a uniformly thick corrugated sheet of ±45° mate- 
rial attached to an equally thick flat sheet of ±45° material with 0° material at specified 
locations in elements 3 and 4. The corrugated configuration is a uniformly thick corru- 
gated sheet of ±45° material with 0° material at specified locations in element 3. In both 
configurations, the thickness of a ±45° sheet is denoted t2(±45°) or simply tg. 


11 




t 3 (0°) 


(a) Hat-stiffened panel. 



(c) Filament orientation pattern. 


Figure 8.- Configurations and design variables. Element identification 
numbers are indicated in figure 8(a). 
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The design variables are the thicknesses t 2 (±45°), tg(0°), and t 4 (0°), and the 

lengths bp b 2 , bg, and b^. In all of the studies, the length L of the panel is 0.76 m 
(30 in.). The upper and lower bounds on the design variables are given in table I, and the 
material properties used for the studies are given in table II. 

Panels Designed for Longitudinal Compression and Shear 

Hat-stiffened panels - longitudinal compression and shear.- A structural efficiency 
diagram for hat- stiffened, graphite -epoxy panels is presented in figure 9. In this diagram 



N 

-X kPa 
L ' 


Figure 9.- Structural eliiciency ol graphite -epoxy, hat-stiffened panels designed 
for pure longitudinal compression and for longitudinal compression and shear. 

the mass index of the panels is shown as a function of the loading index N X ^L for 

two loading cases: (1) pure longitudinal compression and (2) combined longitudinal com- 

pression and shear with N X y^N x = 0.5. The solid curves are for panels designed using 
the present procedure, and the dashed curves are for panels designed using the simplified 
procedure of reference 11. All design variables vary freely and continuously within the 
upper and lower bounds given in table I. For these cases, the present procedure provides 
designs that are up to 14 percent lighter than designs obtained with the simplified pro- 
cedure presented in reference 11. The panels designed for N 
40 percent heavier than panels designed for N v only. 

A 

Corrugated panels — longitudinal compression and shear.- A structural efficiency 
diagram for corrugated, graphite -epoxy panels is presented in figure 10. The curve 


/N x =0.5 


are about 
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Figure 10.- Structural efficiency of graphite -epoxy corrugated panels designed for pure 
longitudinal compression and for longitudinal compression and shear. 
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designations and loadings are similar to those used in figure 9. An important difference 
between the results for the hat-stiffened panels and the results for the corrugated panels 
is that for the corrugated panels the dashed curves fall below the solid curves, indicating 
that the panels designed with the simplified procedure are lighter than the panels designed 
with the present procedure. However, unlike the hat-stiffened panels, calculations with 
VIPASA indicate that corrugated panels designed with the simplified procedure do not 
carry the design load. The adequacy of the simplified analysis is discussed in greater 
detail in appendix B. The panels designed for N 
heavier than panels designed for N x only. 

Hat-stiffened panels - longitudinal compression with stiffener spacing requirement. - 
A common design requirement for stiffened panels is minimum stiffener spacing. The 
effect of imposing such a requirement is shown in figure 11. The curves are for hat- 
stiffened, graphite-epoxy panels supporting a longitudinal compressive loading only. In 
the upper two curves the distance d between stiffeners is equal to 12.7 cm (5.0 in.). 

In the lower two curves the distance d is allowed to vary freely and is always less 
than 12.7 cm. The design trends are the same as those presented in reference 11. The 
mass penalty associated with increasing the stiffener spacing is 30 to 75 percent, 
depending upon the loading. For the hat-stiffened panels with a spacing requirement of 
d = 12.7 cm the present procedure provides designs that are about 7 percent lighter than 
the corresponding designs obtained with the simplified procedure of reference 11 . 


'N = 0.5 are about 60 percent 

X 
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Figure 11.- Structural efficiency of graphite -epoxy, hat-stiffened panels designed 
for longitudinal compression; with and without constraint on distance d 
between stiffeners. 

Effect of Transverse Compressive Load Ny 

The primary inplane loading on a uniaxially stiffened panel is usually a combination 

of N y and N v „, where N Y is in the direction of the stiffeners. The structure con- 
a xy x 

taining the panel is normally designed so that other structural members carry the major 
portion of the transverse load Ny. Small transverse loads in a panel cannot, however, 
be completely eliminated. Studies were carried out with the present procedure to 
examine the effect of small transverse compressive loads Ny on the load-carrying 
ability and on the mass of hat-stiffened panels. 

In the first study, buckling loads were calculated for a panel designed only for a 
longitudinal load N x but subjected to a small transverse load Ny in addition to the 
N x load. The design load is N^L = 689 kPa (100 lbf/in^); Ny is 3 percent of N x . 

The thickness t 2 is 0.56 mm (0.022 in.). Buckling response diagrams for the panel 
with and without the transverse loading are given in figure 12. For this panel, the results 
indicate that with a combined loading of N x and Ny in which N y is only 3 percent 
of N x , the critical value of N x is reduced by 28 percent from its value when Ny = 0. 

It is apparent that a small transverse compressive loading Ny can cause a large reduc- 
tion in the longitudinal buckling load N v . 




Figure 12.- Ratio of buckling load to design load as a function of buckling half -wavelength 
for one hat-stiffened panel designed to support a loading of N^L = 689 kPa 
(100 lbf/in^). In the upper buckling response diagram, the loading is the design 
loading. In the lower diagram, a small transverse loading N y is added to the design 
longitudinal loading N x . t2 = 0.56 mm (0.022 in.). 



Figure 13.- Structural efficiency of graphite-epoxy, hat-stiffened panels designed for 
pure longitudinal compression and for longitudinal and transverse compression. 
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A structural efficiency diagram which shows the increase in mass required to sup- 
port a transverse loading is presented in figure 13. In the lower curve N y is equal to 
zero, and in the upper curve N y is equal to 10 percent of N x . In both curves all 
dimensions, including t 2 > vary freely and continuously. For the range of the loading 
index N x /L considered, the panels with N y y/N x = 0.1 are about 70 percent heavier 
than the panels with no transverse load. It is likely that this increase in mass would be 
diminished if 90° plys were included in the panel. 

Tailoring the Buckling Response of a Panel 

In the buckling response diagram presented in figure 5 the buckling load has a rela- 
tive minimum at X = L/28, which corresponds to a local buckling mode, and another low 
point at X = L which is an overall buckling mode. Other panels in this report have a 
similar buckling response (one or two short-wavelength buckle modes and one or two 
longer wavelength buckle modes are active simultaneously) whether designed by the simpli- 
fied procedure of reference 11 or the present procedure. When additional design require- 
ments such as material strength or extensional and shear stiffness begin to influence the 
design, the number of active buckling modes may decrease. 

In some cases a designer may wish to tailor the buckling response of a panel. For 
example, he may want only one buckling mode to be active at the design load, or he may 
want to put a substantial margin of safety on the buckling load for a given wavelength. 
Because of its generality, the present procedure can be used to tailor the buckling 
response of a panel. 

In order to obtain some insight into this tailoring capability, a heavily loaded hat- 
stiffened panel and a more lightly loaded hat-stiffened panel were considered. The design 
requirements for these two examples are as follows: 


Hat-stiffened 

V L 

t 

2 


1 

Comments 

panel 

kPa 

lbf/in 2 

mm 

in. 

cm 

in. 

Heavily loaded example 

5520 

800 

1.12 

0.044 

Free 

Free 

20-percent margin on buckling for X ^ L/4 

Lightly loaded example 

689 

100 

Free 

Free 

12.7 

j 

5.0 

50-percent margin on buckling for X s L/3 


The buckling response for the heavily loaded example designed with a 20 -percent 
margin on local buckling (X ^ L/4) and the buckling response for the corresponding panel 
designed without the margin on local buckling are presented in figure 14. The panel 
designed with no margin on local buckling (the circular symbols) is buckling critical at 
four wavelengths: X = L, L/2, L/3, and L/20. The 20-percent increase in the local 

buckling load is achieved with a 2.4-percent increase in the mass index. There are no 
stiffener spacing requirements in this example. 
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Figure 14.- Ratio of buckling load to design load as a function of buckling half -wavelength 
for two hat-stiffened panels designed for N x y/L = 5520 kPa (800 lbf/in 2 ); t 2 = 1.12 mm 
(0.044 in.). 
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Figure 15.- Ratio of buckling load to design load as a function of buckling half -wave length 
for two hat-stiffened panels designed for N X/ /L = 689 kPa (100 lbf/in 2 ); d = 12.7 cm 
(5.0 in.). 
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Figure 16.- Mass index as a function of required local buckling load N l /L and required 
overall buckling load Nq/L for graphite -epoxy, hat-stiffened panels with and without 
design requirements on distance d between stiffeners. Required buckling loads 
given in kPa (lbf/in^). 


The buckling response for the more lightly loaded example designed with a 
50-percent margin on overall buckling (X § L/3) and the buckling response for the corre- 
sponding panel designed without the margin on overall buckling are presented in fig- 
ure 15. Because of the stiffener spacing requirement, the local buckling wavelength for 
each of these panels is larger than for panels without that additional requirement. The 
50-percent increase in the overall buckling load is achieved with a 3.2-percent increase 
in the mass index. 


Additional studies were carried out with the more lightly loaded panel - both with 
and without a requirement on the stiffener spacing d. In these studies the mass penalty 
associated with increasing the local and/or overall buckling loads is examined. The 
results are presented in the form of carpet plots in figure 16. The stiffener spacing d 
is equal to 12.7 cm (5.0 in.) in the upper plot and varies freely in the lower plot. The 
quantities that are varied in the carpet plots are the required local buckling load 

X g L/4) and the required overall buckling load ^Nq^L, X ^ L/3^ . The vertical 


axis is the mass index 


W/A 

L 


The thickness t2 varies freely. 


All of these tailoring studies indicate that if a stiffened panel has no stiffener- 
spacing requirement, the local buckling load can be increased with a relatively small 
mass penalty. If, on the other hand, a stiffened panel has a stiffener-spacing requirement 
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that has a strong effect on the design, as it does in this case (see fig. 11), then- the overall 
buckling load can be increased with a relatively small mass penalty. 

Although imperfections are not considered in this report, it is known that imperfec- 
tions can cause a substantial reduction in the buckling load of a panel. In many cases the 
reduction is caused by an interaction between imperfections and various possible buckling 
modes. Wavelengths of both imperfections and buckling loads are important. If a 
designer knows the type of imperfection that will occur during fabrication and/or the type 
of damage that is likely to occur during service, he may be able to tailor the buckling 
response of a panel to reduce the detrimental effects of the imperfections. 

CONCLUDING REMARKS 

A procedure is described for designing uniaxially stiffened panels made of composite 
material and subjected to combined inplane loads. The procedure uses a rigorous buck- 
ling analysis and nonlinear mathematical programing techniques. Computation efficiency 
is attained by using approximate analysis techniques based on Taylor series expansions. 
The procedure is very general and is applicable to panels with arbitrary cross sections. 
The procedure has been applied in a computer program for designing hat- stiffened and 
corrugated panels. 

Design studies, all of which consider graphite -epoxy panels, are also presented. In 
one set of studies, generalized mass-strength charts are developed for hat-stiffened and 
corrugated panels subjected to combined longitudinal compression and shear. For the 
hat-stiffened panels, the procedure provides designs that are up to 14 percent lighter than 
those obtained with a previously published procedure based on a simplified buckling analy- 
sis. For the corrugated panels, the studies show that buckling modes involving torsional 
rolling may cause the simplified procedure to be unconservative. Buckling modes 
involving torsional rolling are treated accurately in the present procedure. 

In other studies, the effects of combined longitudinal and transverse compression 
are investigated for hat-stiffened panels. Small transverse compressive loads severely 
reduce the load-carrying capacity of longitudinally stiffened panels. For transverse 
loads about one -tenth of longitudinal loads, panel mass must be increased by as much as 
70 percent above the mass of panels with no transverse load. 

Finally, the capability to tailor the buckling response of hat-stiffened panels under 
uniaxial compression is explored. Results indicate that relatively small mass increases 
could be associated with putting large margins on either (a) local buckling if no stiffener - 
spacing requirement is imposed, or (b) overall buckling if a minimum stiffener-spacing 
requirement is imposed. Carpet plots are presented for a limited loading range to 
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suggest mass penalties associated with buckling-response tailoring for panels with and 
without stiffener -spacing constraints. 

The present study demonstrates that a panel design procedure with a high-quality 
buckling analysis and with complete generality of constraints is practical. Such pro- 
cedures can be used to avoid failure from complex buckling modes and to determine 
mass and proportions of panels for multiple design load conditions and constraints. 

Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
March 24, 1977 
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THE VIPASA ANALYSIS CODE 

This appendix gives a brief description of the VIPASA analysis code in terms of the 
analysis -design procedure and examples presented herein. Two main topics are con- 
sidered: (1) a general discussion of VIPASA with emphasis on structural modeling, and 

(2) conservatism in the VIPASA buckling analysis in the case of a shear loading and/or 
anisotropy. A more complete description of VIPASA is given in references 3, 4, and 12. 

General Discussion of VIPASA 

The computer program VIPASA calculates the buckling loads of structures com- 
prised of flat rectangular plate elements connected together along their longitudinal edges. 
The response of each element is obtained using an exact solution of the thin-plate equa- 
tions with anisotropic terms included. The analysis connects the individual plate ele- 
ments and maintains continuity of the buckle pattern across the intersection of neighboring 
plate elements. Individual plate elements may be isotropic, orthotropic, or anisotropic. 
The only limitations are that the plates must be uniform along their length and that the 
plate laminates must not exhibit bending-stretching coupling or stretching-shearing 
coupling. (Bending-twisting coupling is allowed.) Laminates that are midplane sym- 
metric and balanced satisfy these requirements. The structure may be subjected to any 
combination of longitudinal, transverse, and shear loads that are constant along the length. 

In the VIPASA analysis, a local coordinate system is defined for each individual 
plate element. An example is shown in figure 17. The X, Y^, Z^ axes define the 
local coordinate system in the local longitudinal, transverse, and lateral directions, 
respectively. The buckling displacements u, v, and w are defined in this local coor- 
dinate system and are also shown in figure 17. The X, Y, Z axes denote the overall -panel 
coordinate system. 

In VIPASA, the buckling boundary conditions along the edges x = 0 and x = L 

(fig. 18) are always simple support - that is, u is unrestrained, v = w = 0, and 

w is unrestrained. The unconnected longitudinal edges may, however, have various 

boundary conditions. In the examples presented in this report, the boundary conditions 

along the unconnected longitudinal edges are either simple support or symmetric. 

(Calculations were made with both boundary conditions.) Symmetric boundary conditions 

along the unconnected longitudinal edges are defined as u is unrestrained, v = 0, w is 

unrestrained, and w =0. The net effect of the individual plate-element boundary con- 

,y fi 

ditions is to produce identical u,v,w boundary conditions in the X,Y, Z, respectively, 
coordinate system associated with the overall panel. 
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Figure 17.- Local coordinate system for panel element and overall 
coordinate system for stiffened panel. 


In the VIPASA model for the hat-stiffened and corrugated panels considered in this 
report, one stiffener section is idealized as an assemblage of flat plates, then eight of 
these substructures are joined together for the final model. A sample cross section is 
shown in figure 18. Because eight stiffeners are used to model all of the panels discussed 



Figure 18.- Eight-stiffener mathematical model of hat-stiffened panel. 


in this report, symmetric and simple support boundary conditions along the unconnected 
longitudinal edges give virtually identical values for the buckling load. The potential for 
buckle widths (y -direction) to extend across more than eight stiffeners increases with the 
presence of shear loadings and that potential is even greater for transverse loadings. In 
the results presented in figures 9 and 13, checks were made at a few design points to 
insure that eight stiffeners were adequate and that the minimum buckling load had been 
found. 
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(b) Mathematical model used in VIPASA. 

Figure 19.- Approximate physical model for hat-stiffened panel and 
corresponding mathematical model used by VIPASA. 


A detailed model of one repeating section of a hat-stiffened panel is shown in fig- 
ure 19(a). Values for plate stiffnesses supplied to VIPASA are based on the fact that 
adjacent elements are connected through nodes at the element center lines. When plates 
do not align center line to center line, which is the case in figure 19(a), eccentricities or 
offsets must be supplied to VIPASA in order for the panel to have the proper stiffnesses. 
The eccentric connections for the hat-stiffened panel are illustrated in figure 19(b). The 
nodes are indicated by dots, and the offset distances EZ1, EZ2, and EZ3 represent 
eccentricities in the z-direction. Element 1 is given a positive offset EZ1 so that the 
top surface of the panel is continuous and flat through its width. At the juncture of ele- 
ment 2 with element 4, element 2 is given a negative eccentricity EZ2. Where element 2 
joins element 3, element 2 is given a positive eccentricity EZ3. The sign convention used 
in figure 19(b) is consistent with that used in VIPASA. 

Conservative Analysis in the Case of Shear Loading and/or Anisotropy 

In order to obtain an exact solution to the panel buckling problem, the following 
buckling mode shape is assumed in the VIPASA analysis for each plate element: 


24 



I 



For orthotropic plates without shear loading, F is real and, hence, node lines are 
straight, perpendicular to the x-direction, and spaced X apart. In this case, simple 
support boundary conditions are achieved at the boundaries x = 0 and x = L. With 
shear and/or anisotropy present, F is complex and node lines are skewed and not 
straight, but they are still spaced X apart. An example is shown in figure 20. If X is 

PANEL BOUNDARIES 

NODAL LINES DETERMINED 

BY V1PASA 



Figure 20.- Nodal lines determined by VIPASA in the case of 
a shear loading and/or anisotropy. 

small compared with panel length, the unrestrained buckle pattern can develop in the cen- 
tral region of the panel. The VIPASA results are, in this case, sufficiently accurate. 

As X approaches L, the VIPASA results can be quite conservative. The reason is as 
follows. As can be seen in figure 20, the skewed nodal lines do not coincide with the panel 
edges x = 0 and x = L. The additional constraint arising from the requirement that the 
nodal lines at x = 0 and x = L coincide with the panel edges results in long -wavelength 
buckling loads that are, in many cases, appreciably higher than those determined by 
VIPASA, where the nodal lines are not required to conform with the edges. Calculations 
have shown that for long-wavelength buckling modes the effect of anisotropy is minimal. 
Anisotropy, therefore, causes negligible conservatism in a VIPASA analysis. The pres- 
ence of a shear loading can, however, lead to very conservative results for X equal to L. 

Because of VIPASA’ s conservatism in the case of long-wavelength buckling if a 
shear load is present, a modified analysis procedure is used for the case X = L. The 

25 



APPENDIX A 


modified analysis is based on an interaction formula used in reference 11. In the case 
of a combined longitudinal and shear loading ^N x and N xy j the loading is considered 
to be critical if 



.1 


(A2) 


The critical longitudinal compressive load N is calculated with VIPASA. The crit- 

x cr 

ical shear load N is calculated using smeared orthotropic theory which is the same 

■ cr 

approach used in reference 11. It is emphasized that this modified analysis procedure is 
used only for the case A = L. 
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TESTING THE ADEQUACY OF THE SIMPLIFIED ANALYSIS-DESIGN 
PROCEDURE OF NASA TN D-8257 

In NASA TN D-8257, two forms of buckling are postulated. One form is an overall 
mode, and the other form is a local mode. Buckling of the panel is assumed to be pre- 
vented if these two forms of buckling are prevented. 

In the overall buckling mode, the panel is assumed to buckle as a wide column. All 
the orthotropic stiffnesses are smeared. This failure mode is often referred to as Euler 
buckling. In the local buckling mode, each element of the panel cross section is treated 
as a narrow orthotropic plate that is simply supported along its lines of attachment to 
adjacent elements. Buckling of each element is treated separately (i.e., one can speak 
of buckling of element 1, buckling of element 2, etc.). This approach to local buckling 
ignores continuity of buckling-mode shape. The advantage of this type of analysis pro- 
cedure is that closed-form solutions exist for the buckling modes that are postulated. 

The disadvantage is that some buckling modes - in particular those involving joint dis- 
placement (rolling modes, for example) - are not considered. The analysis approach 
used in reference 11 is a traditional approach that has been used in many other panel 
optimization procedures. 

In contrast with the simplified analysis of reference 11, the VIPASA analysis code 
provides a high-quality buckling analysis that considers all buckling modes and insures 
continuity of the buckling pattern across the intersection of neighboring plate elements. 
Anisotropy in the element bending stiffnesses can also be taken into account. The major 
known shortcoming of the VIPASA code is the skewed boundary node lines (rather than 
straight boundary node lines) determined by VIPASA in the case of a shear loading and/or 
anisotropy. That characteristic and a modification to the analysis are discussed in 
appendix A. 

In this appendix the adequacy of the analysis -design procedure of reference 11, 
including the orthotropic assumption, is checked using the present procedure. The 
checking is done with two levels of analysis. First, designs obtained with the procedure 
of reference 11 are analyzed with VIPASA under the assumption of orthotropy. In these 
cases differences in buckling loads cannot be ascribed to anisotropic effects. Second, 
the same designs are analyzed with VIPASA using the anisotropic stiffnesses Djg and 
Dgg which couple bending behavior with torsional behavior. (See, for example, eq. (AT) 
in ref. 11.) Finally, the masses of the panels designed using the procedure of reference 11 
are compared with the masses of panels designed using the present procedure. 
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The analysis and design comparisons that are presented in table HI fall under two 
main headings: "Analysis Comparisons" and "Design Comparisons." Each of these 

headings is subdivided into two subheadings depending upon the analysis used in VIPASA. 

In the analysis comparisons, the data presented are the ratios of the buckling load calcu- 
lated by VIPASA to the buckling load calculated using the procedure of reference 11. The 
panels being analyzed were designed using the procedure of reference 11; therefore, the 
buckling load calculated by the procedure of reference 11 is the design load. The buckling 
half -wavelength X determined by VIPASA is also presented. In the design comparisons, 
the data presented are the ratios of the mass index of panels designed using the present 
procedure to the mass index of panels designed using the procedure of reference 11. 

Hat-Stiffened Panels 

Longitudinal compression N x .- Five hat-stiffened panels that were designed using 

the procedure of reference 11 to support longitudinal compressive loads N x were 
analyzed with VIPASA. Panels were also designed using the present procedure so that 
panel masses could be compared. The results are summarized in table III. 

More valuable information on the buckling characteristics, including the anisotropic 
effect, can be obtained by examining the buckling response over a range of buckling wave- 
lengths. A buckling response diagram obtained with the VIPASA code for a hat-stiffened 
panel having t£ = 0.56 mm (0.022 in.) and designed with the procedure of reference 11 
for a loading of N x /L = 689 kPa (100 lbf/in2) is presented in figure 21. The circular 

"ST TT 



Figure 21.- Ratio of buckling load to design load as a function of buckling half -wavelength 
for a hat-stiffened panel designed using the procedure of reference 11 for a loading of 
of N x //L = 689 kPa (100 lbf/in^); tg = 0.56 mm (0.022 in.). 


28 



APPENDIX B 


symbols indicate the calculated buckling loads for a given half -wavelength obtained using 
an orthotropic VTPASA analysis. The square symbols are for a VIPASA analysis including 
anisotropic D 16 and D 26 terms. The critical overall buckling load occurs at X = L, 
and the critical local buckling load occurs at X = L/16. The buckling mode shapes for 
these two values of X are similar to those presented in figures 6 and 7, respectively. 

The main difference is that in the panel of figure 21 the local buckling mode does not 
involve deformation of the hat cap (element 3 in fig. 8). 

Three conclusions can be drawn from the results presented in figure 21 and from 
the first set of studies in table in. First, for the cases examined, anisotropy has its 
main effect on buckling modes with wavelengths less than the critical wavelengths for 
local buckling. Therefore, for these cases, anisotropy has little impact on the buckling 
load for longitudinal compressive loadings. Second, the simplified analysis of refer- 
ence 11 appears to be adequate* for hat-stiffened panels loaded by longitudinal compres- 
sion. However, it would be inadequate if the wavelengths for local buckling and strong 
anisotropic effects coincide. (Calculations indicate that the D^g and D 2 g terms for 
elements 1 and 2 of fig. 8 provide almost all of the anisotropic effects for the short- 
wavelength buckle modes.) Third, for hat-stiffened panels with a longitudinal compres- 
sive loading N x , mass savings on the order of 2 to 9 percent are possible with the present 
design procedure. 

Longitudinal compression and shear N x and N x y- Three hat- stiffened panels that 

were designed using the procedure of reference 11 to support longitudinal compressive 

and shear loads N x and N X y were analyzed with VIPASA. Panels were also designed 

using the present procedure so that panel masses could be compared. The results are 

summarized in table III. The positive and negative values of N„./n indicate the sign 

xy/ a 

of the shear loading N x y. Because the analysis of reference 11 is based on orthotropic 
theory, it cannot distinguish between positive and negative values of N . In contrast, 
since VIPASA can take into account the anisotropic D^g and D 20 terms VIPASA can 

2 

distinguish between positive and negative values of N X y. In some cases presented in 
table III, the anisotropic effects are large. Panels designed with an orthotropic analysis 
can have an anisotropic buckling load substantially less for a negative value of N xy than 
for a positive value of N xy . Also, the mass of a panel designed for a negative value 

*As used here and elsewhere in this appendix, the word "adequate" means that the 
panels designed using the procedure of reference 11 do not buckle when analyzed using the 
VIPASA code. In many cases the more sophisticated analysis -design procedure pre- 
sented herein provides designs that are substantially lighter than designs obtained with 
the procedure of reference 11. 

^The effect of a sign change for is the same as the effect of changing the sign 

of the D^g and D 20 terms. Positive values of D 16 and D 20 are used in all the 
results presented. 
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of Nxy can be substantially greater than the mass of a panel designed for a positive 
value of Njjy. 

A buckling response diagram provides a clearer picture of the buckling character- 
istics including local buckling and anisotropic effects. The buckling response diagrams 
for the panel designed for N x yi. = 689 kPa (100 lbf/in 2 ), ^ xy /^ x = 0.1 and for the 
panel designed for N X/ /l = 207 kPa (30 lbf/in 2 ), ^ xy ^ x = are presented in fig- 
ures 22 and 23, respectively. Both panels were designed with the procedure of refer- 
ence 11. The circular symbols are for a VIPASA analysis that assumes complete 
orthotropy. The square and triangular symbols are for VIPASA analyses that account for 
D lg and D 26 terms. The square symbols indicate positive values of N xy and the 
triangular symbols indicate negative values. The solid symbol at A = L is for the mod- 
ified analysis described in appendix A. (See eq. (A2).) 

Considering first the orthotropic VIPASA analysis indicated by the circular sym- 
bols, the analysis of reference 11 is adequate for the short -wavelength buckles. However, 
when the VIPASA analysis includes the D lg and D 2g terms there is a substantial dif- 
ference in the short-wavelength buckling loads depending on the sign of the shear loading. 
These effects cannot be detected with the analysis of reference 11. 

It can be concluded that for combined longitudinal compression and shear loadings 
the local buckling analysis in reference 11 is adequate for hat-stiffened panels made up of 
orthotropic elements. When anisotropic effects are large, as they can be in a four -ply 
laminate, the analysis of reference 11 is inadequate, and a modified local buckling analysis 
that accounts for anisotropic effects must be used. 

An analysis similar to that of reference 11 which takes into account the effects of 
the D^ g and D 2g terms in the local buckling analysis for elements 1 and 2 is easily 
developed. The technique can be used if the thickness t 2 is fixed during the synthesis. 
With t 2 fixed, the D.j matrix shared by elements 1 and 2 is fixed. The local buckling 
loads of elements 1 and 2 can, therefore, be expressed as 


f N x 'i = ~/ — ~To 

(Bl) 

\ x crjk (b k ) 2 

/N \ = K 9 

(B2) 

\ ^cr/k ( b k ) 2 


where b^ is the width of element k, and C and K are constants, each of which can 
be determined by carrying out a single anisotropic buckling analysis with a code such as 
VIPASA. Sample calculations indicate that the resulting analysis produces conservative 
results for hat-stiffened panels. 
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Figure 22.- Ratio of buckling load to design load as a function of buckling half -wavelength 
for a hat-stiffened panel designed using the procedure of reference 11 for a loading 
of N x /L = 689 kPa (100 lbf/in 2 ); N xy /N x = 0.1; t 2 = 0.56 mm (0.022 in.). 





L 

Figure 23.- Ratio of buckling load to design load as a function of buckling half -wavelength 
for a hat-stiffened panel designed using the procedure of reference 11 for a loading 
of N x /L = 207 kPa (30 lbf/in 2 ); N xy/ /N x = 1.0; t 2 = 0.56 mm (0.022 in.). 
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Corrugated Panels 

Longitudinal compression N x> - Five corrugated panels that were designed using 

the procedure of reference 11 to support longitudinal compressive loads N x were 
analyzed using VTPASA. Panels were also designed using the present procedure so that 
panel masses could be compared. The results are presented in table in and figures 24 
and 25. 




Figure 24.- Ratio of buckling load to design load as a function of buckling half -wavelength 
for a corrugated panel designed using the procedure of reference 11 for a loading of 
N x /L = 1380 kPa (200 lbf/in^); tg = 0.56 mm (0.022 in.). 



Figure 25.- Buckling-mode shape for panel whose buckling response diagram 

is shown in figure 24. 
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It can be seen in table ID that deficiencies in the simplified analysis of reference 11 
exist for this configuration even when comparisons are made between orthotropic analyses. 
The analysis of reference 11 appears to be adequate in the lightly loaded region of a given 
thickness t 2 , but is very unconservative in the heavily loaded region of a given thickness. 
(These regions can be seen in fig. 19 of ref. 11.) The buckling response diagram for the 
panel designed for N x /L = 1380 kPa (200 lbf/in 2 ) with t 2 -0.56 mm (0.022 in.) is pre- 
sented in figure 24. The buckling mode shape for the X = L buckling mode of that panel 
is presented in figure 25. Although the simplified (ref. 11) analysis for compression 
buckling of corrugated panels appears to be adequate for short -wavelength buckling, it is 
not adequate for the longer wavelength buckling X = L and L/2 where, according to the 
buckling-mode shape presented in figure 25, the corrugated stiffeners are rolling or 
twisting. Twisting modes cannot be predicted by the analysis of reference 11. As in the 
hat-stiffened panel, anisotropic effects are small in the case of corrugated panels loaded 
only in longitudinal compression. 

Longitudinal compression and shear N x and N X y.- Two corrugated panels 

designed using the procedure of reference 11 to support longitudinal compressive and 
shear loads N x and N X y were analyzed using VIPASA. Panels were also designed 
using the present procedure so that panel masses could be compared. The results are 
summarized in table III and figure 26. 



Figure 26.- Ratio of buckling load to design load as a function of buckling half -wavelength 
for a corrugated panel designed using the procedure of reference 11 for a loading of 
N X/ /L = 689 kPa (100 lbf/in 2 ); N xy/ /N x = 0.1; t 2 = 0.56 mm (0.022 in.). 
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The buckling response diagram obtained with the VIPASA code for the t 2 = 0.56 mm 
(0.022 in.) panel designed with the reference 11 analysis for a loading of N x /L = 689 kPa 
(100 lbf/in2) ; N /n = 0.1 is presented in figure 26. The figure indicates that the 
simplified analysis would be adequate for this corrugated panel if it were made up of 
orthotropic elements. As in the case of the hat-stiffened panel, anisotropy affects mainly 
the short-wavelength buckling loads. 
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TABLE I.- LIMITS ON DESIGN VARIABLES 


Dimension 

tg(0°) ? minimum 
t^(0°) ? minimum 
t^O 0 ), maximum 
t 4 (0°), maximum 
b^ minimum 
b2 ? minimum 
minimum 
b^ 7 minimum 


Hat-stiffened 


panel 


cm 

in. 

0 

0 

0 

0 

0.318 

0.125 

.635 

.250 

2.03 

.80 

1.02 

.40 

2.03 

.80 

1.02 

.40 


Corrugated 


panel 


cm 

in. 

0 

0 

0 

0 

0.318 

1 0.125 

Not applicable 

2.03 

.80 

1.02 

.40 

2.03 

.80 

Not applicable 



TABLE H.- PROPERTIES OF GRAPHITE -EPOXY MATERIAL 


USED IN SAMPLE CALCULATIONS 


[ten. denotes tension; comp, denotes compressio^ 


Symbol 

Value in 
SI Units 

Value in U.S. 
Customary Units 

Density and elastic properties 

p 

1520 kg/m^ 

0.055 lbm/in^ 

E i 

145 GPa 

21 x 10 6 psi 

E 2 

16.5 GPa 

2.39 x 10 6 psi 

G 12 

4.48 GPa 

0.65 x 10 6 psi 

m 12 

0.314 

0.314 

^21 

0.037 

0.037 


Allowable stresses 

cr^ (ten.) 

1240 MPa 

180 ksi 

cr^ (comp.) 

-1240 MPa 

-180 ksi 

cr^ (ten.) 

55.2 MPa 

8 ksi 

&2 (comp.) 

-207 MPa 

-30 ksi 

a 

t 12 

82.7 MPa 

12 ksi 


Allowable strains 

(ten.) 

0.0087 

0.0087 

(comp.) 

-0.0087 

-0.0087 

e| (ten.) 

0.00475 

0.00475 

€2 (comp.) 

-0.01764 

-0.01764 

a 

*12 

0.01846 

0.01846 



TABLE HI.- RATIOS OF BUCKLING LOADS AND MASS INDICES FOR GRAPHITE -EPOXY PANELS ANALYZED 
AND DESIGNED BY PROCEDURE OF REFERENCE 11 AND BY PRESENT PROCEDURE 







| Analysis comparisons 1 

Design comparisons 






1 N \ 


( W/ A \ 


Design 

loading 


to (±4 5°) 

\ cr / V1PASA 

\ L / present procedure 


£ 


1 IK \ 


/ W/A \ 







\ cr / reference 11 

\ L / reference 11 

N x 

L 

N xy 



VIPASA 

VIPASA 

Present 

procedure 

orthotropic 

Present 

procedure 

anisotropic 

kPa 

lbf 

in2 

N x 



orthotropic 

anisotropic 





Hat-stiffened panels, N x 



345 

50 

0 

0.056 

0.022 

1.04 X = L 

1.03 X = L/14 

0.97 

0.98 

689 

100 

0 

.056 

.022 

1.03 X = L/16 

1.02 X = L/16 

.95 

.97 

1380 

200 

0 

.056 

.022 

1.01 A = L 

1.01 A = L 

.91 

.96 

2070 

300 

0 

.056 

.022 

.98 X = L/2 

.98 X = L/2 

.96 

.98 

2070 

300 

0 

.112 

.044 

1.00 X = L 


.97 






Hat- stiffened panels, N and N 

x XV 



689 

100 

+0.1 

0.056 

0.022 

1.04 X = L/16 

1.05 X = L/16 

0.93 

0.93 



-.1 




.89 X = L/38 


1.00 

1380 

200 

+.1 

.056 

.022 

1.02 X = L/12 

1.03 X = L/12 

.94 

.92 



-.1 




1.01 X = L/12 


.98 

207 

30 

+1.0 

.056 

.022 

1.03 X = L 

1.03 X = L 

.85 

.83 

— 


-1.0 




.84 X = L/43 


1.46 





Corrugated panels, N x 



345 

50 

0 

0.056 

0.022 

1.03 X = L 

1.00 X = L/19 

0.96 

0.98 

689 

100 

0 

.056 

.022 

.96 X = L/2 

.96 X = L/2 

1.00 

1.01 

1380 

200 

0 

.056 

.022 

.66 A = L 

.66 A = L 

.97 

.98 

2760 

400 

0 

.056 

.022 

.57 A = L 

.57 A = L 

1.11 

1.13 

1380 

200 

0 

.112 

.044 

.82 X = L/2 


1.06 






Corrugated panels, N x 

and N xy 



345 

50 

+0.1 

0.056 

0.022 

1.03 X = L ' 

1.03 X = L 

0.96 

0.92 



-.1 




.90 X = L/29 


1.02 

689 

100 

+.1 

.056 

.022 

1.00 X = L 

1.00 X = L 

.86 

.84 



-.1 




.97 X = L/38 


1.01 


* Designs obtained with procedure of reference 11. 
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